In silico discovery of epitopes of gag and env proteins for the development of a multi-epitope vaccine candidate against Maedi Visna Virus using reverse vaccinology approach.

Maedi Visna Virus (MVV) causes a chronic viral disease in sheep. Since there is no specific therapeutic drug that targets MVV, development of a vaccine against the MVV is inevitable. This study aimed to analyze the gag and env proteins as vaccine candidate proteins and to identify epitopes in these proteins. In addition, it was aimed to construct a multi-epitope vaccine candidate. According to the obtained results, the gag protein was detected to be more conserved and had a higher antigenicity value. Also, the number of alpha helix in the secondary structure was higher and transmembrane helices were not detected. Although many B cell and MHC-I/II epitopes were predicted, only 19 of them were detected to have the properties of antigenic, non-allergenic, non-toxic, soluble, and non-hemolytic. Of these epitopes, five were remarkable due to having the highest antigenicity value. However, the final multi-epitope vaccine was constructed with 19 epitopes. A strong affinity was shown between the final multi-epitope vaccine and TLR-2/4. In conclusion, the gag protein was a better antigen. However, both proteins had epitopes with high antigenicity value. Also, the final multi-epitope vaccine construct had a potential to be used as a peptide vaccine due to its immuno-informatics results.

Ovine PAPPA2 gene coding variants are linked to decreased fecal egg shedding in native Turkish sheep naturally infected with gastrointestinal nematodes.

In this study, the association between PAPPA2 coding variants and gastrointestinal (GI) nematode fecal egg count (FEC) score in adult Turkish sheep was investigated. For this purpose, the FEC score was determined in adult sheep from six breeds: Karacabey Merino (n = 137), Kivircik (n = 116), Cine capari (n = 109), Karakacan (n = 102), Imroz (n = 73), and Chios (n = 50). Sheep were classified as shedders or non-shedders within breeds and flocks. The first group was the fecal egg shedders (> 50 per gram of feces), and the second group was the no fecal egg shedders (≤ 50 per gram of feces). The exon 1, exon 2, exon 5, exon 7, and a part of 5'UTR of the ovine PAPPA2 gene were genotyped by Sanger sequencing of these two groups. Fourteen synonymous and three non-synonymous single-nucleotide polymorphisms (SNPs) were found. The non-synonymous SNPs, D109N, D391H, and L409R variants, are reported for the first time. Two haplotype blocks were constructed on exon 2 and exon 7. The specific haplotype, C391G424G449T473C515A542 on the exon 2 that carries the 391H variant, was tested against four other common haplotypes. Our results indicate that C391G424G449T473C515A542 haplotype was significantly associated with fecal egg shedding status in adult Turkish sheep (p-value, 0.044).

Excessive replacement changes drive evolution of global sheep prion protein (PRNP) sequences.

Sheep prion protein (PRNP) is the major host genetic factor responsible for susceptibility to scrapie. We aimed to understand the evolutionary history of sheep PRNP, and primarily focused on breeds from Turkey and Ethiopia, representing genome-wise ancient sheep populations. Population molecular genetic analyses are extended to European, South Asian, and East Asian populations, and for the first time to scrapie associated haplotypes. 1178 PRNP coding region nucleotide sequences were analyzed. High levels of nucleotide diversity driven by extensive low-frequency replacement changes are observed in all populations. Interspecific analyses were conducted using mouflon and domestic goat as outgroup species. Despite an abundance of silent and replacement changes, lack of silent or replacement fixations was observed. All scrapie-associated haplotype analyses from all populations also showed extensive low-frequency replacement changes. Neutrality tests did not indicate positive (directional), balancing or strong negative selection or population contraction for any of the haplotypes in any population. A simple negative selection history driven by prion disease susceptibility is not supported by the population and haplotype based analyses. Molecular function, biological process enrichment, and protein-protein interaction analyses suggested functioning of PRNP protein in multiple pathways, and possible other functional constraint selections. In conclusion, a complex selection history favoring excessive replacement changes together with weak purifying selection possibly driven by frequency-dependent selection is driving PRNP sequence evolution. Our results is not unique only to the Turkish and Ethiopian samples, but can be generalized to global sheep populations.

A novel 2 bp deletion variant in Ovine-DRB1 gene is associated with increased Visna/maedi susceptibility in Turkish sheep.

Visna/maedi (VM) is a multisystemic lentivirus infection of sheep that affecting sheep industry across the globe. TMEM154 gene has been identified to be a major VM-associated host gene, nevertheless, a recent study showed that the frequency of the VM-resistant TMEM154 haplotypes was very low or absent in indigenous sheep. Thus, the present study was designed to determine other possible co-receptors associated with VM. For this purpose, DRB1 gene, which is renowned for its role in host immune response against various diseases was targeted. A total number of 151 case-control matched pairs were constructed from 2266 serologically tested sheep. A broad range of DRB1 haplotype diversity was detected by sequence-based genotyping. Moreover, a novel 2 bp deletion (del) in the DRB1 intron 1 was identified. For the final statistic, the sheep carrying VM-resistant TMEM154 diplotypes were removed and a McNemar's test with a matched pairs experimental design was conducted. Consequently, it was identified for the first time that the 2 bp del variant is a genetic risk factor for VM (p value 0.002; chi-square 8.31; odds ratio 2.9; statistical power 0.90) in the dominant model. Thus, negative selection for 2 bp del variant could decrease VM infection risk in Turkish sheep.

Association of TLR2 haplotypes encoding Q650 with reduced susceptibility to ovine Johne's disease in Turkish sheep.

Ovine Johne's disease (OJD) is caused by Mycobacterium avium subsp. paratuberculosis (MAP) and carries a potential zoonotic risk for humans. Selective breeding strategies for reduced OJD susceptibility would be welcome tools in disease eradication efforts, if available. The Toll-like receptor 2 gene (TLR2) plays an important signaling role in immune response to MAP, and missense variants are associated with mycobacterial infections in mammals. Our aim was to identify and evaluate ovine TLR2 missense variants for association with OJD in Turkish sheep. Eleven TLR2 missense variants and 17 haplotype configurations were identified in genomic sequences of 221 sheep from 61 globally-distributed breeds. The five most frequent haplotypes were tested for OJD association in 102 matched pairs of infected and uninfected ewes identified in 2257 Turkish sheep. Ewes with one or two copies of TLR2 haplotypes encoding glutamine (Q) at position 650 (Q650) in the Tir domain were 6.6-fold more likely to be uninfected compared to ewes with arginine (R650) at that position (CI95 = 2.6 to 16.9, p-value = 3.7 × 10-6). The protective TLR2 Q650 allele was present in at least 25% of breeds tested and thus may facilitate selective breeding for sheep with reduced susceptibility to OJD.

Evaluation of CD109, PCP4 and SEMA3D genes for their association with Ovine Johne's disease in Turkish sheep.

Johne's disease is a chronic, contagious, zoonotic disease that affects numerous species including livestock and sometimes humans. The disease is globally distributed in sheep populations and caused by Mycobacterium avium Subsp. paratuberculosis (MAP). A previous genome-wide association study identified single nucleotide polymorphism (SNP) markers associated with OJD serostatus in CD109, PCP4, and SEMA3D genes. Our aim was to evaluate the same markers for association with OJD seroprevalence in Turkish sheep in a retrospective matched case-control study. The serological status for OJD in 1801 sheep was determined for four native and four composite breeds from three research flocks. One hundred eleven matched case-control pairs were constructed according to breed type and age from 1750 comingled ewes reared in the same environment. A Single Nucleotide Primer Extension (SNuPE) assay was designed to genotype PCP4-Intron 1, PCP4-3'UTR, SEMA3D, CD109-intron 2 and CD109-intron 8 markers and a McNemar's test was performed on the matched pairs. An association with these five markers was not detected with the OJD serostatus in Turkish sheep (power of detection, 0.95; odds ratio >3; McNemar's p < .05). Thus, a wider search may be needed to identify any major underlying genetic risk factors for OJD in Turkish sheep.

Documentation of extensive genetic diversity in the Ovar-DRB1 gene in native Turkish sheep.

Indigenous breeds have a high level of genetic diversity that might contribute to develop animal breeds with desired traits such as disease resistance and high productivity. Major histocompatibility complex (MHC) is a key component of adaptive immune system and consists of highly polymorphic genes that take part in adaptive immune response and disease resistance. Exploring and understanding the effect of polymorphisms in MHC could be beneficial to future animal breeding strategies. In this study, we sequenced the highly polymorphic Exon2 of the ovine DRB1 gene using Sanger sequencing to explore the diversity of this gene in six indigenous Turkish sheep breeds and two crossbreeds. In total, 894 haplotypes from 447 sheep were investigated, and 69 different haplotypes including 27 novel ones were identified. Among the identified haplotypes there were common and breed specific haplotypes. There was a relatively high diversity of the alleles within indigenous breeds. Allelic diversity patterns were mostly associated with geographical differences. The results of this study highlight the genetic variation within indigenous breeds which has important implications for biodiversity and the adaptability of breeds to specific environments. There is value to further studies which include other genomic regions and traits, and these could guide breeding strategies.

Association between ovine Toll-like receptor 4 (TLR4) gene coding variants and presence of Eimeria spp. in naturally infected adult Turkish native sheep.

Coccidiosis caused by Eimeria spp. is a protozoan disease prevalent in farm animals, and it is responsible for serious economic losses especially in young animals. It has been popular to breed disease-resistant animals due to the concern about food safety, animal welfare, and public health. Toll-like receptor (TLR) gene family plays a key role in the innate immune system participating in host-antigen interaction, therefore, they are candidate genes for breeding disease-resistant animals. In the present study, possible genetic associations between TLR4 gene coding variants and the presence of Eimeria spp. in adult Turkish sheep were investigated. For this purpose, the presence of Eimeria spp. in fecal samples from six native Turkish sheep were determined, and approximately 1450 bp region in the 3rd exon of the ovine TLR4 gene was sequenced. Ten nonsynonymous and four synonymous single nucleotide polymorphisms (SNPs) were detected in the targeted region. Statistical analyses revealed that the SNP at the codon at 356th position encoding Leucine instead of Phenylalanine (F356L) was significantly associated with the presence of Eimeria spp. It was found that the individuals carrying at least one Leucine amino acid sequence at this position have 2.3-fold more risk for the presence of Eimeria spp.

Detecting fecal egg count (FEC) for gastrointestinal nematodes of adult Turkish sheep with different scrapie related PRNP haplotypes.

Scrapie is a transmissible spongiform encephalopathy caused by prions and leads to neurodegeneration in the Central Nervous System (CNS) of sheep and goats. Genetic resistance/susceptibility to scrapie is well studied and it is known that the variations of 136th, 154th and 171st codons at the ovine PRNP gene have a major effect on the development of the disease. Many studies demonstrated that selection for PRNP genotypes has not influenced other performance traits, nevertheless, there is a knowledge gap about the possible link between the PRNP gene and the status of the other important diseases that affect the sheep population worldwide. In the present study, we tested whether there is an association between scrapie-related PRNP genotypes and fecal egg count (FEC) of gastrointestinal nematodes in seven adult Turkish sheep breeds. For this purpose, FEC scores of studied sheep (n = 253) were determined and the same animals were genotyped for the PRNP gene. Finally, an association analysis was performed for scrapie resistant (ARR), susceptible (VRQ), and wild-type (ARQ) haplotypes. Based on our statistical analysis, it is concluded that PRNP genotypes have no positive or negative effect on the FEC scores of adult sheep.

Five novel PRNP gene polymorphisms and their potential effect on Scrapie susceptibility in three native Ethiopian sheep breeds.

BACKGROUND: Classical scrapie susceptibility in sheep has been linked to three polymorphisms at codon 136, 154, and 171 in the prion protein gene (PRNP) whereas atypical scrapie susceptibility is related to polymorphisms at codon 141. Many other variants over the length of the PRNP have been reported. Some of the variants may play crucial roles in fighting against the emergence of a new form of scrapie disease. Scrapie surveillance, scrapie associated genotyping and PRNP characterization studies have been conducted across the globe. However, such in-depth studies have never addressed the African continent's sheep breeds. Therefore, genotyping native Ethiopian sheep breed's PRNP gene has socioeconomic and scientific merits. This study aimed to identify PRNP variants in three native Ethiopian sheep breeds and their potential effect on scrapie susceptibility. RESULTS: Five novel variants were identified in the PRNP gene of three native Ethiopian sheep breeds. Four non-synonymous heterozygous substitutions i.e. H99Q (CAC-- > CAA), H99L (CAC-- > CTA), A116E (GCA-- > GAA), A116T (GCA-- > ACA), and one synonymous N103 N (AAC-- > AAT) were detected. In addition to the novel variants, polymorphisms at codon 126,127,138,142,146,231, and 237 were also identified. The haplotype ARR was observed in Menz and Afar breeds at frequencies of 0.02 and 0.05 respectively. Neither ARR/ARR nor VRQ/VRQ genotypes were identified in the population under study. CONCLUSION: Two of the novel variants at codon 99 and 103 that are placed closer to the proteinase K cleavage site and the variant at codon 116 in the palindrome region along with variants at codon 127 in glycine repeat domain may influence the conformational flexibility of prion protein. The rarity of ARR haplotype and the abundance of 141 L variant demonstrated that the present study population was less resistant to classical scrapie and less predisposed to genotype associated atypical scrapie. This study provides a valuable dataset that can be potentially integrated into selective breeding strategies during interbreeding, crossbreeding and help to take precautionary measures against scrapie.

Novel Variations in Native Ethiopian Goat breeds PRNP Gene and Their Potential Effect on Prion Protein Stability.

Scrapie is a lethal neurodegenerative disease of sheep and goats caused by the misfolding of the prion protein. Variants such as M142, D145, S146, H154, Q211, and K222 were experimentally found to increase resistance or extend scrapie incubation period in goats. We aimed to identify polymorphisms in the Afar and Arsi-Bale goat breeds of Ethiopia and computationally assess the effect of variants on prion protein stability. In the present study, four non-synonymous novel polymorphisms G67S, W68R, G69D, and R159H in the first octapeptide repeat and the highly conserved C-terminus globular domain of goat PrP were detected. The resistant genotype, S146, was detected in >50% of the present population. The current study population showed a genetic diversity in Ethiopian goat breeds. In the insilico analysis, the R68 variant was predicted to increase stability while S67, D69, and H159 decrease the stability of prion protein. The new variants in the octapeptide repeat motif were predicted to decrease amyloidogenicity but H159 increased the hotspot sequence amyloidogenic propensity. These novel variants could be the source of conformational flexibility that may trigger the gain or loss of function by prion protein. Further experimental study is required to depict the actual effects of variants on prion protein stability.

Nucleotide and octapeptide-repeat variations of the prion protein coding gene (PRNP) in Anatolian, Murrah, and crossbred water buffaloes.

Resistance to bovine spongiform encephalopathy (BSE) that is significantly associated with insertion/deletion (indel) polymorphisms at two loci (putative promoter and intron 1) on the prion protein gene (PRNP) in cattle has been well documented. Studies suggest that the insertion alleles are related to BSE resistance. Until recently, BSE has never been reported in water buffaloes (unlike cattle). Previous studies have demonstrated that the PRNP gene in water buffalo consists mostly of insertion alleles at both loci; nevertheless, whether or not water buffaloes are genetically resistant to BSE and the role of indel polymorphisms in their resistance status is not clear. We examined the coding region of PRNP to determine the nucleotide and octapeptide-repeat (octarepeats) variations of Anatolian, Murrah and Murrah × Anatolian (M × A) water buffaloes. Three synonymous single nucleotide polymorphisms (SNP) at positions 126, 234, and 285, and a non-synonymous SNP at position 322 (G108S) were detected. Triplet G/A/T base substitutions were observed at position 126 and two additional genotypes, T/A and T/G, at this position were determined. We also found six octarepeats that indicated the presence of the wild-type PRNP6 allele in the coding region. To the best of our knowledge, this is the first report of the T/A and T/G genotypes in water buffaloes.

Investigation of the prion protein gene (PRNP) polymorphisms in Anatolian, Murrah, and crossbred water buffaloes (Bubalus bubalis).

Bovine spongiform encephalopathy (BSE) of the cattle is the outstanding disease among other transmissible spongiform encephalopathy (TSEs). It can be transmitted from the cattle to a human and causes a new variant of the Creutzfeldt-Jakob disease (CJD). It is known that prion protein coding gene (PRNP) plays a major role in the TSE susceptibility or resistance in some species. Recent researches demonstrated that the insertion (in) and deletion (del) polymorphisms within promoter and intron 1 region of the PRNP related to BSE susceptibility in cattle. In contrast to cattle, BSE has never been reported in water buffalo; hence, PRNP polymorphisms may be an explanation for buffalo resistance to BSE. The aim of this study was to evaluate allele, genotype, and haplotype frequencies of the PRNP promoter and intron 1 insertion/deletion (indel) polymorphism in healthy Anatolian, Murrah, and Murrah × Anatolian crossbred buffaloes. According to our findings, there were no deletion alleles at two mentioned loci. All studied buffaloes were monomorphic and have carried in/in haplotypes which are considered as the most resistant genotype to BSE.